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Abstract 



An effective theory approach to heavy quarkonia decays based on the 1/rriQ 
expansion is introduced. Its application to decays in which the two heavy quarks 
' annihilate is discussed. 
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§ : 1 Introduction 

Heavy Quark Effective Theory (HQET) [|I|] has turned out to be a very successful approach 
to describe systems with a single heavy quark. It is based on the infinite mass limit of 
Q-i! QCD which serves as a starting point to perform a systematic expansion in Kqcd/^q 
and a s {mo) using the methods of effective field theory. 

Presently we are at the threshold of the discovery of flavoured "doubly heavy" systems, 
i.e. states consisting of two heavy quarks (such as the _B c 's and also baryons with b = c = 
±1, b = ±2 or c = ±2), while quarkonia like systems (the f/>'s and the T's) are in the 
meantime quite well known. 

Motivated by these expectations the question arises whether one may set up a similar 
effective theory approach for systems with two (or even more) heavy quarks based on the 
I/tuq expansion of QCD. It turns out that one can not use the static limit for two heavy 
quarks if their velocities differ only by an amount of order 1/uiq, i.e. vv'—l ~ Kqcd/^q- 
The static limit breaks down and one is forced to include at least the kinetic energy into 
the leading order dynamics; in other words one has to use a non-relativistic approximation 
instead of the static limit. Such approximations have been formulated in the QED context 
(NRQED) some time ago; Two slightly different approaches (NRQCD § and HQQET §) 
have recently been studied in QCD. In this talk the general features of such an effective 
theory is outlined and its application to this class of decays is considered in which the 
two heavy quarks annihilate. 



2 Structure of HQQET/NRQCD 



In order to set up a heavy quarkonium effective theory one starts from the I/tuq expansion 
of the QCD Lagrangian and the corresponding expansion of the fields 
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The first terms of these two expansions define the static limit, which has been successfully 
applied to systems with a single heavy quark. In order to describe a system with more 
than one heavy (anti) quark one has to write down the same expansion ([]]) for each heavy 
quark. However, in the static limit for a state with two or more heavy quarks one runs 
into problems with diverging phases and "complex anomalous dimensions", which are 
considered in detail in 0. 

In order to cure this problem one has to choose the unperturbed system such that these 
phases are already generated by the leading order dynamics, i.e. instead of the static limit 
one has to use the non-relativistic Lagrangian. For a system with a heavy quark and a 
heavy antiquark one then starts from 

C ^h^\ivD)h^-hi-\ivD)h^ + K u ^ 1 = ^+)|^!/ i (+) + ^-)|^!/ i H (2) 

2m Q 2m Q 

where we have assumed for simplicity that the two quarks have the same mass; the case 
of unequal mass is obvious. 

Most of the success of HQET is due to heavy quark flavour and spin symmetry. 
However, once one uses (0) the symmetries are somewhat different for HQQET. First of 
all, (0) depends on the mass through the kinetic energy term; consequently the states 
will depend on itlq in a non-perturbative way and heavy flavour symmetry is lost. On 
the other hand, (|2|) does not depend on the spins of the two heavy quarks so there is a 
spin symmetry which is larger than in HQET because we have two heavy quark spins; 
the resulting symmetry is an SU(2) <g> SU(2) corresponding to separate rotations of the 
two spins. 

For the case of heavy quarkonia all states fall into spin symmetry quartets which 
should be degenerate in the non-relativistic limit. In spectroscopic notation 25+1 lj these 
quartets consist of the states 

[nH, n 3 4-i n 3 £ e n 3 £ i+1 ] (3) 

For the ground states the spin symmetry quartet consists of the t/q (the 0~ state) and 
the three polarization directions of the Tq (the 1~ state). 

The heavy quarkonia spin symmetry restricts the non-perturbative input to a calcu- 
lation of processes involving heavy quarkonia. Of particular interest are decays in which 
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the heavy quarks inside the heavy quarkonium annihilate. The annihilation is a short dis- 
tance process that can be calculated perturbatively in terms of quarks and gluons, while 
the long distance contribution is encoded in certain matrix elements of quark operators. 
Logarithmic dependences on the heavy quark mass may be calculated by employing the 
usual renormalization group machinery. 



3 Annihilation Decays of Heavy Quarkonia 

The starting point to calculate processes like tjq — > light hadrons, tjq — > 7 + light hadrons, 
or the corresponding decays of the Tq states is the transition operator T for two heavy 
quarks which annihilate into light degrees of freedom. This will in general be bilinear in 
the heavy quark fields, such that 

T(X, = (-i)Q(X + £)K(X, $Q(X - (4) 

where K(X,£) involves only light degrees of freedom and X and £ correspond to the cms 
and relative coordinate respectively. We identify the field Q with the quark and Q with 
the antiquark, so we shall make the large scale itlq explicit by redefining the fields as 

Q(x) = exp(-zm Q (fx))(5[ +) (x), Q(x) = exp(-im Q (vx))Q { ~\x) (5) 

This corresponds to the usual splitting of the heavy quark momentum into a large part 
ttiqv and a residual piece k. Inserting this into (f|) this yields 

T=(-i)exp[-i2m Q vX)Qi-\X + 0K(X,t)Qi + \X-0 (6) 

The inclusive decay rate for the decay of a quarkonium $ — > light degrees of freedom is 
then given by 

T = (^\\[ d A Xd 4 £d 4 CT(X, £)T\0, + h.c.l |tf ) (7) 



The next step is to perform an Operator Product Expansion (OPE) for the non-local 
product of the quark field operators. This expansion will yield four-quark operators of 
increasing dimension starting with dim-6 operators. The increasing dimension of these 
operators will be compensated by inverse powers of the heavy quark mass, so generically 
the rate takes the form 



n,i \ m Qj 



where n — 6, 7, • - • is the dimension of the operator and i labels different operators with the 
same dimension. The coefficients C{0 < f'\ii) are related to the short distance annihilation 
process and hence may be calculated in perturbation theory in terms of quarks and gluons. 
Once QCD radiative corrections are included, the C(C>l: n \ fi) acquire a dependence on the 
renormalization scale /i which is governed by the renormalization group of the effective 
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theory. The rate T is independent of \i and hence the \i dependence of C{0\ , //) has to 
be compensated by a corresponding dependence of the matrix elements. 

The non-perturbative contributions are encoded in the matrix elements of the local 
four-quark operators, and the mass dependence of these operators is also expanded in 
powers of I/ttlq and thus the remaining itlq dependence of the matrix elements is only 
due to the states. In terms of the uiq independent static fields 

Q(+)( x ) = h M (x) + 0(l/m Q ) Qi~\x) = h(~\x) + 0(l/m Q ) 
one has in total four dim-6 operators 

A f) = j i (+)^ Ch (-) j l {-) l5Ch {+) and A f) = lW lji ch { - ) U^^Ch^ 

where C is a color matrix, where one has the two possibilities C®C=l®loiC®C = 
rpa rpa^ These operators do not mix under renormalization, all anomalous dimensions 
vanish. 

There are no dim-7 operators, since these are all proportional to (ivD) and can be 
rewritten in terms of dim-8 operators by the equations of motion. At dim-8 one finds 30 
local operators 

{iD,)U + h,Ch^\ [{iD»)U-hsCh^\ 

(iD»)U + ^ 5 Ch^} [h^ l5 C{t B M )/i (+) ] + h.c. 
(iD^h^^Ch^} [#->C(i#)/i (+) ] + h.c. 
(i^)A (+ VC/i H ] [h { ~h»C(i D)xh {+) ] + h.c. 
(iD x )U +) ^Ch^} [U'^xCii 5)^ (+) ] + h.c. 

h {+ h,C(i D),h^\ [h^ l5 C(i B)^ (+) ] 
h^C{i$)h^\ [h^C(ifi)h^\ 
h {+ hxC(i 5)„/i (_) ] [h { -^ x C(i B)^ (+) ] 
h i+) l^C{i D)xh { - ] } [U~^ x C(i DYh {+) ] 

U + h b Ch^} [U-^ 5 C(i Dfh {+) ] + h.c. 
h^-fCh^] [h^C(i]p){i D)^h (+) ] + h.c. 
fe (+ VC/i H ] [h^Cii D)^ip)h (+) ] + h.c. 
hM-fCh^] lh { -h»C(i D) 2 h {+) ] + h.c. 
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In addition to these contributions one has also non-local terms originating from single 
insertions of the Lagrangian of order 1/ttlq and from double insertions of the Lagrangian 
of order I/tuq. Under renormalization the local dim-8 operators do not mix; only the 
double insertion of the kinetic energy operator of order I/uiq mixes into some of the 
above operators. Denoting this contribution as 



TP = i ~Y 1 J dVyTlAPK^K^y)} (9) 

one obtains in one-loop renormalization group improved perturbation theory two sets of 
equations for the coefficients of the operators with the spin structure 75 <8> 75 

C{E?\y) = C{E?\m Q ) - 33_ 8 2n/ g ( T i (1) ^Q) ln ^ 
C{Bf\^ = C{Bf\m Q ) - QQ 24 C(T?\m Q )lnr, 
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<7(£>?\a0 = C{Df\m Q ) - QO C(Ti 1 \m Q )\nr ] + f 33 ^ 2n C(Ti 1 \m Q )\n V 

C(e[ s \^) = C(E[ s \m Q ) - y 3 7 3 -^C( T1 i 1) ^Q) ln ^ 

where r] = (a s (jj)/a s (mQ)). Furthermore, the coefficients C{T^ c \mQ) are the same as 
the ones for the dim-6 operators C(A < f : \ mq) since the kinetic energy operator is not 
renormalized. 

The second set of equations is for the operators with spin structure 7 M <S> 7 M and due to 
heavy quarkonia spin symmetry one obtains the same equations; all other renormalization 
group equations are trivial. 

A calculation of an annihilation decay then involves to calculate the C(Oi n \fi) at 
the scale fi = rriQ by matching the effective theory to full QCD. Once this is done, one 
may run down to some small scale fi of the order of the "binding energy" of the heavy 
quarkonium, thereby resumming the well known logarithms of the form ln(mg//z) that 
appear in the calculations of decay rates of heavy p-wave quarkonia. As an example, in 
0] the decay t\q — > 7 + light hadrons is studied in HQQET. 

The matrix elements of these operators are non-perturbative quantities, which are 
constrained by heavy quarkonia spin symmetry. In order to exploit this symmetry, one 
may use the usual representation matrices for the spin singlet and spin triplet quarkonia 



H^v) = VMP+75 for S = , H 3 (v) = VMP + e for S = 1 (10) 

where M ~ 2mq is the mass of the heavy quarkonium and P + = (1 + f)/2. Using this 
one finds for the matrix elements of the dim-6 operators 

(tf|fcWrC7iHfcHr'C7iW|tf) = a {c \n,£)G, with G = Tt(H 2s+1 T)Tt(T'H 2s+1 ) (11) 
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Thus for each n and £ and for each color combination one finds a single parameter for 
both the spin singlet and spin triplet quarkonium. 
Correspondingly one finds for the dim-8 operators 

(*\[hMrc(i B M )/i ( - } ] (iD u )[h^rch^} + h. c .|*) 

(V\[h™rC(i D„)(i D v )h^\ [h^T'Ch^] + h.c.\^) 

For fixed values of n and £ one finds that eight parameters are needed to describe the 
matrix elements of the dim-8 operators. 

These matrix elements are non-perturbative, but from vacuum insertion one suspects 

a (1) (n, 0) ~ | ^(0) | 2 > a (1) (n, t) for £ ^ , a (1) (n, 0) > a (8) (n, for all n, £ 
d (1) (n, 1) ~ Ki(0) | 2 > d {1 \n,£) for £ ^ 1 , 1) > d {8) {n,£) for all n,£ 

e (1) (n, 0) ~ Re K o (0)i?; o (0)] > e (1) (n, for £ ^ , e (1) (n, 0) > e (8) (n, £) for all n, £ 

where R n i(r) is the radial wave function of the quarkonium. The same reasoning yields 
the expectation that (n, £) and (n, £) are small compared to the coefficients that 
are non-vanishing in vacuum insertion. 
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